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Theoretical Characterization of the Potential Energy Surface 

for NH + NO 
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ELORET Institute 
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Abstract. The potential energy surface (PES) for NH + NO has been characterized 
using complete active space self-consistent field (CASSCF) gradient calculations to 
determine the stationary point geometries and frequencies followed by CASSCF/ 
internally contracted configuration interaction (CCI) calculations to refine the ener- 
getics. The present results are in qualitative accord with the BAC-MP 4 calculations 
of Melius and Binkley, but there are differences as large as 8 kcal/mol in the detailed 
energetics. Addition of NH to NO on a 2 A' surface, which correlates with N 2 + OH 
or H + N 2 0 products, involves barriers of 3.2 kcal/mol (trans) and 6.3 kcal/mol 
(cis). Experimental evidence for these barriers is found in the work of Bohmer et al. 

The 2 A" surface has no barrier to addition, but does not correlate with products. 
Surface crossings between the barrierless 2 A" surface and the 2 A' surface may be 
important. Production of N 2 OH products is predicted to occur via a planar 
saddle point of 2 A! symmetry. This is in accord with the preferential formation 
of H(A') A doublet levels of OH in the experiments of Patel-Misra and Dagdigian. 
Addition of NH 1 A to NO is found to occur on an excited state surface and is 
predicted to lead to N 2 0 product as observed by Yamasaki et al. 

“Mailing Address: NASA Ames Research Center, Moffett Field, CA 94035. 
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I. Introduction 


The reaction of NH with NO is important in the combustion of nitrogen containing 
fuels and in the thermal De-NO* process in which NO is converted to N 2 by the 
addition of NH 3 [1-4]. The latter process is important in controlling NO* emissions 
from jet engines. 

Considering the X 3 E _ ground state of NH, there are two exothermic product 
channels. 


NH X 3 E _ + NO —> N 2 + OH (1) 

NH X +NO ^ H + N 2 0 (2) 

Experimental estimates of the product branching ratio vary widely. Yamasaki et al. 
[5] concluded that reaction (1) accounted for 100 % of the products, while Mertens 
et al. [6] reported that reaction (1) accounted for only 19 ± 10 % of the products. 

The PES for reaction (1) [7,8,9] involves an intermediate HNNO complex which 
undergoes a 1,3-hydrogen shift with a barrier to form an NNOH species which is 
unstable and dissociates to N 2 + OH. The HNNO complex can also dissociate over 
a barrier to H + NNO. The reverse process 


H + NNO -+ N 2 + OH 

( 3 ) 

H + NNO -» NH + NO 

( 4 ) 


has also been studied. Reaction (3) can occur via a direct process in which H 
attaches to the O end of N 2 0, or via an indirect process in which the H attaches to 
the end N of N 2 0 and then undergoes a 1,3-hydrogen shift leading to products via 
the same saddle point as for reaction (1) [8]. 
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The interned energy distribution of the OH product for reaction (1) has been 
studied by Patel-Misra and Dagdigian [10] in a crossed beam experiment. The 
OH v=l to v=0 vibrational population ratio was found to be 0.30 ± 0.06 and the 
average rotational energy was found to be 25 ± 1 kJ/mol (5.97 ±0.23 keal/mol). 
These authors note that this is a rather small amount of internal excitation given 
that reaction (1) is exothermic by 408 kJ/mol (97.5 keal/mol). Attempts to model 
the internal energy distribution in the OH product by phase space theory lead 
to more excitation than observed, suggesting that the HNNO complex may not 
be long enough lived to exhibit statistical behaviour. In this study the internal 
energy distribution in the N 2 was not measured and it was presumed that most 
of the excess energy went into relative translational energy of the products. More 
recently, Bohmer et al. [11] have studied the reaction of N 2 O with hot H atoms 
formed by photolysis of an N 2 O-HI complex. This study is in general agreement 
with the work of Patel-Misra and Dagdigian with respect to the OH internal energy 
distribution, but the OH Doppler shift profiles in these experiments require a large 
amount of internal vibrational excitation of the N 2 ( ~ 25,000 cm -1 (71.5 keal/mol) 
for OH v=0). This result is explained by a model in which the HNNO complex has 
an elongated NN bond and the 1,3-hydrogen shift is assumed to occur rapidly with 
respect to motion of the heavy atoms leaving the N 2 product vibrationally excited. 
A more detailed understanding of the product energy distributions and branching 
ratios in this system requires more detailed information about the PES, which is 
reported in this paper. 

Unlike many other reactions of NH, where the 1 A state is more reactive, the NH 
+ NO reaction is found to have comparable rates for the ground X 3 E - state and the 
low lying 1 A state of NH [5]. Yamasaki et al. [5] find a strong product specificity 
depending on the NH spin state. The ground state of NH leads to N 2 + OH as the 
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predominant product, while the 1 A state of NH also leads to H + N 2 O. In addition, 
NO is able to quench NH 3 A to NH 3 E - . 

NH 1 A + NO -» H + N 2 0 (5) 

NH 1 A + NO—*NH X 3 E~ + NO (6) 

In the work of Yamasaki et al. the experimental results are analyzed in terms of a 
model based on a calculated PES by Melius and Binkley [7] and Fueno et al. [9]. 
In this picture NH (*A) + NO correlates with a 2 A' surface which also correlates 
with H + N 2 O, while NH (X 3 S~) + NO correlates with a HNNO complex on the 
2 A" surface. As discussed below, this model is substantially oversimplified. In fact 
NO + NH in the 3 E“ and 1 A state lead to six surfaces which are strongly coupled 
and the reaction of the 1 A state of NH is a carbene insertion process which occurs 
on an excited state surface. Preliminary calculations on the excited state surfaces 
axe also reported here. 

There have been several previous studies of the PES for NH + NO. The most 
extensive of these is by Melius and Binkley [7] and used the BAC-MP4 method 
with a 6-31G** basis set. More recently Fueno et al.[9] have also carried out multi 
reference configuration interaction (MRCI) studies of this system using smaller basis 
sets. The present calculations are similar in spirit to these earlier studies, but use 
larger basis sets and more extensive configuration interaction. Thus, the present 
results are expected to be more reliable than the earlier studies. Qualitative features 
of the potential energy surfaces are discussed in Sec. II, the computational method 
is discussed in Sec. Ill, the results are presented in Sec. IV, and the conclusions 
are given in Sec. V. 

II. Qualitative Features. 
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Combining the 3 E“ ground state of NH with the 2 II ground state of NO and 
keeping a plane of symmetry leads to 2,4 A' and 2 ’ 4 A" surfaces. (It should be noted 
here that for four atoms the reaction may occur without symmetry; however, the key 
stationary points on the surface are found to have Cj symmetry and this symmetry 
is assumed in the following discussion.) 

For the 2 A" surface, the singly occupied NO 27 t orbital is in the plane of the 
supermolecule and addition of NH to NO occurs with no barrier. 




However, a doubly occupied 0 2p orbital is in the plane of the supermolecule and 
1,3-hydrogen migration is unfavorable. Thus, approach on this surface would not 
be expected to lead to N 2 + OH products. Also this surface can not lead to H + 
N 2 0, since this product channel is of 2 A' symmetry. 

For the 2 A' surface, the singly occupied NO 2ir orbital is perpendicular to the 
plane of the supermolecule and the NH is adding to an NO tt bond. Thus, a barrier 
to addition is expected on this surface. However, a singly occupied O 2p orbital is in 
the plane of the supermolecule and 1,3-hydrogen migration leading to NNOH may 
occur. This process is expected to have a barrier since an NH bond is being broken 
as the OH bond is formed. The NNOH species is not bound and the 1,3-hydrogen 
migration leads directly to N 2 + OH. 

Fig. 1 shows a schematic of the potential energy surface for NH + NO. Focusing 
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first on the features discussed above, it is seen that NH + NO may form HNNO 
in the cis or trans orientation with no barrier on the 2 A" surface, but there are 
barriers on the 2 A' surface. HNNO corresponds to a stable planar minimum on the 
PES and the 2 A' state is lower than the 2 A" state. This leads to entrance channel 
surface crossings between these two surfaces, which may lead to surface hopping 
between the 2 A" surface, which forms with no barrier, and the 2 A' surface, which 
has an entrance channel barrier but leads to products. cis-HNNO on the 2 A' surface 
may undergo a 1,3-hydrogen shift leading to N 2 + OH products. 

From Fig. 1 it is also seen that H + NNO correlates with the 2 A' surface. 
Addition of H to the 0 end of N 2 O leads directly to N 2 + OH wdth a large barrier. 
An alternative indirect pathway involves addition of H to the end N of N 2 O which 
leads to cis-HNNO with a smellier barrier than addition to the 0 end. This pathway 
then connects to N 2 + OH by the same 1,3-hydrogen shift saddle point involved 
in the NH + NO reaction. The indirect pathway has a lower overall barrier and 
should be the predominate pathway at lower temperatures. 

Addition of the low-lying 1 A excited state of NH to NO involves a geometry in 
which the N atom of NH approaches near the midpoint of the NO bond with the 
NH bond approximately perpendicular to the plane formed by the N atom of NH 
and the NO bond. Combining the 3 E — and *A states of NH with the 2 H state of 
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NO leads to six potential energy surfaces. Of these the lowest two are the 2 A 7 and 
2 A” surfaces which correlate with the 3 E~ ground state of NH. Of the remaining 
four surfaces three are repulsive and one is attractive. The electronic configuration 
of the attractive state is: 



(where the doubly occupied N 2p orbital perpendicular to the plane represents the 
NH bond pair of NH) This configuration is analagous to the configuration involved 
in insertion of 0( 1 D) into H 2 [12]. This configuration can also correlate with H + 
N 2 0, since dissociation of the NH bond leads to a configuration with four a 7 and 
four a" (with respect to the NNO plane) electrons in i r orbitals on the N 2 0. It is 
probable that this excited surface accounts for the production of N 2 0 observed by 
Yamasaki et al. [5] in the reaction of NH 3 A with NO. 

III. Computational Details. 

Two different basis sets were used in this work. For the CASSCF gradient calcu- 
lations the polarized double zeta set of Dunning and Hay [13] was used. The basis 
set for N and 0 is a (9s5p)/[3s2p] basis augmented by a single set of 3d functions 
with exponents of 0.80 and 0.85 for N and 0, respectively. The H basis is (4s)/ [2s] 
augmented with a single set of 2p functions with exponent 1.00. The basis set used 
in the Cl calculations is the Dunning correlation consistent triple zeta double po- 
larization basis set [14]. This basis is [4s3p2dlf] for N and 0 and [3s2pld] for H 
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and is described in detail in Ref. 14. 

The CASSCF calculations had 13 electrons distributed among 8 a' and 3 a" 
orbitals. The active electrons correspond to the electrons depicted in Eqns. 7 and 
8 plus the 2s electron pair on the central N. The remaining 10 electrons are the two 
pairs of N Is electrons, the 0 Is and 2s electron pairs, and the N 2s electron pair on 
the end N. One N 2s pair has to be active to describe certain regions of the surface 
(e.g. H + N 2 O) where the N 2s electrons on the central N are used in bonding. 

All but the N Is and 0 Is electrons are correlated in the subsequent CCI calcu- 
lations. A selected reference list was used with a coefficient cutoff of 0.04. 

The CASSCF/gradient calculations used the SIRIUS/ABACUS system of pro- 
grams [15], while the CCI calculations were carried out with MOLPRO [16,17]. 
Most of the calculations were carried out on the NASA Ames Cray Y-MP. Some of 
the CCI calculations were carried out on the NAS facility YMP. 

IV. Discussion. 

Before discussing the results for the HN + NO reaction, it is useful to consider re- 
sults for the molecular species involved. Table I compares geometries and harmonic 
vibrational frequencies from CASSCF/gradient calculations for NH, NO, N 2 O, OH, 
and N 2 with experiment. The results shown in Table I axe typical of results obtained 
at stable minimum points on a PES using CASSCF/gradient methods. Here it is 
seen that the bond lengths axe too long by « 0.02 A but the vibrational frequen- 
cies are in fair agreement with experiment. This result contrasts with the results 
of Hartree-Fock theory, which typically gives bond lengths which are somewhat 
shorter than experiment and frequencies which axe larger than experiment. It has 
been common practice to scale SCF frequencies by ~ 0.9 to account for this defect 
in the SCF results. The CASSCF frequencies on the other hand are approximately 
correct and are used here without scaling. 
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Tables Ha and lib show comparisons of CASSCF/gradient and externally con- 
tracted Cl [18] calculations for the HN 2 and H0 2 surfaces. Once again at the 
minima on the PES’s, corresponding to HN 2 and H0 2 , the geometries and frequen- 
cies are in fair agreement with the more accurate CCI results. ( The bond lengths 
from the CASSCF calculations are slightly too long as noted above.) At the saddle 
points the CASSCF results are poorer. In both cases the saddle point geometry 
is tighter in the CASSCF calculations than in the CCI calculations, as evidenced 
by a shorter HN/HO bond length, a larger bending frequency, and a larger imagi- 
nary frequency at the saddle point. It is also seen from Table II that the binding 
energies are smaller in the CASSCF calculations. These results are typical, based 
on experience with other systems. For HN 2 it has been shown that addition of 
a correction term, which is a slowly varying function of the NH bond length, to 
the CASSCF surface results in a reasonable approximation to the CCI surface [19]. 
This result suggests that the CASSCF surface is mainly distorted along the reaction 
coordinate, but is a better representation in the directions orthogonal to the reac- 
tion coordinate. Thus, one way to correct the CASSCF surface near a saddle point 
would be to compute Cl energies along the minimum energy path on the CASSCF 
surface. This approach has actually been carried out for the reaction of J CH 2 + 
H 2 0 [20], which is a barrierless process. However, for the present system the saddle 
points involve barriers and are expected to be less sensitive to small changes of the 
geometry along the reaction coordinate. Thus, the approach used here is to carry 
out CCI calculations at the stationary point geometries obtained at the CASSCF 
level. 

Tables Ilia and Illb give computed frequencies and geometries for the 2 A' and 
2 A" states of HNNO, respectively. Table IV gives the computed CCI energies at 
each of the stationary points, the zero-point energies, and the relative energies ( 
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including zero-point energy). The same energetics are shown graphically in Fig. 1. 

From Table Ilia it is seen that addition of NH to NO on the 2 A' surface occurs 
with loose saddle points i.e. the bond lengths and frequencies are close to those 
for NH and NO. The saddle points for adding H to the end N and 0 of N 2 O show 
somewhat more distortion of bond lengths, ryvjv elongated by 0.03 A for H-NNO 
and tno elongated by 0.07 A for NNO-H. 

Cis and trans HNNO on the 2 A ( surface have similar bond lengths. Both the NN 
and NO bond lengths of HNNO on the 2 A' surface are shorter than on the 2 A" 
surface. The NN bond length is shorter by 0.10 A and 0.12 A for cis and trans 
HNNO, respectively. This effect results from the greater NN 7 r bonding on the 
2 A' surface which is evident from Eqn’s (7) and (8). The shorter NO bond on the 
2 A' surface is probably due to smaller non-bonded repulsions between the doubly 
occupied O 2p orbital and the NN bond orbitals. 

The saddle point for 1,3-hydrogen migration in cis-HNNO leading to N 2 + OH 
product (denoted NNOH in Table Ilia), by contrast with the other saddle points 
discussed above, would be described as a tight saddle point. Here the NN bond 
is elongated by 0.11 A (compared to a comparable N 2 calculation, Table I). As 
discussed by Bohmer et al. [11] (see section I), this elongation of the NN bond is 
believed to be responsible for the production of vibrationally excited N 2 product. 

Table V compares the geometries and frequencies obtained in the present calcu- 
lations with the BAC-MP4 results of Melius and Binkley [7,21]. It should be noted 
that the frequencies reported by Melius and Binkley have been scaled downward 
by 12 %. For cis-HNNO the CAS bond lengths are somewhat longer than those 
obtained with the BAC-MP4 method. This is expected from the discussion above. 
The frequencies are in fair accord but there appears to be a tendency for the CAS 
frequencies to be larger than the scaled UHF frequencies. For the saddle points the 
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CASSCF geometries are again in reasonable accord with the IJHF geometries (tak- 
ing into account the tendency for CAS to slightly overestimate bond lengths). One 
interesting difference is that CAS leads to a larger imaginary frequency for NNO-H 
compared to NNOH, while UHF shows the opposite trend. This could have some 
significant dynamical consequences, because in Ref. 8 a one dimensional tunneling 
correction was made based on the imaginary frequency. 

Table VI compares the CAS/CCI energetics with those obtained using the BAC- 
MP4 method. Here there is sporadic agreement. In particular, the NNO-H and 
NNOH saddle point energies agree to within 2 kcal/mol with the present calcula- 
tions, but the H-NNO saddle point is 6.7 kcal/mol lower in the present calcula- 
tions. A peculiar result is that the BAC-MP4 calculations find cis-HNNO below 
trans-HNNO which is the opposite of what is found in the present calculations. It 
should be noted here that the BAC-MP4 method makes rather large corrections 
10 kcal/mol for single bonds and « 20 kcal/mol for double bonds) to the en- 
ergies. The present calculations use no empirical corrections yet yield results of at 
least as good accuracy. The incorrect ordering for cis and trans HNNO is puzzling 
because the correct ordering is obtained with CASSCF and CCI. Still the present 
results tend to support the utility of the BAC-MP4 method, though the present 
large scale Cl calculations are undoubtedly more reliable. 

We now turn to a more detailed discussion of the computed potential energy 
surface as shown in Fig. 1 and based on the energetics given in Table IV. As shown 
in Fig. 1 there are barriers to the addition of NH to NO on the 2 A' surface. The 
computed barrier heights are 6.3 kcal/mol for cis approach and 3.2 kcal/mol for 
trans approach. These barriers are expected based on the discussion in Section II 
and are similar in origin to barriers for addition of H to NO [22]. Experimental 
verification of these barriers is obtained in experiments of Bohmer et al. [11] which 
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measure the threshhold for formation of NH 3 S~ in the reaction of hot H atoms 
with NO formed by photolysis of an HI-NO complex. These experiments indicate a 
threshhold occuring for a photolysis wavelength between 250 and 260 nm whereas 
the energy of the lowest vibrational level of NH 3 E~ + NO occurs at « 270 nm, a 
difference of « 4-9 kcal/mol. It should be noted here that H + NO correlates with 
the 2 A' surface and dissociation to NH + NO therefore should involve a barrier. 

Addition on the 2 A" surface involves no barrier but this surface does not correlate 
with H + N 2 0 nor, from the discussion in Section II, is it likely to lead to OH -f- N 2 
products. Because the 2 A' surface drops more rapidly than the 2 A" surface, surface 
crossings are expected in the NH + NO entrance channel region. These surface 
crossings may allow addition to occur on the barrierless 2 A” surface with subsequent 
surface crossing to the 2 A' surface which connects to products. This process could 
reconcile the apparent lack of a barrier to reaction (1) with the expectation that the 
products pass through the planar NNOH saddle point, which has 2 A' symmetry. 
Evidence that the products arise via a saddle point of 2 A' symmetry comes from 
the experiments of Petel-Misra and Dagdigian [10] in which a marked preference for 
H(A') A doublet levels is observed in the OH product. 

The reverse reaction of H + N 2 O i.e reactions (3) and (4) has been studied by 
Marshall et al. [8] and by Bohmer et al. [11]. The present results are consistent 
with the discussion in Ref. 8. The production of N 2 + OH can proceed via a direct 
process in which the H adds to the O end of N 2 O. The computed barrier height for 
the direct pathway is 18.0 kcal/mol. The indirect pathway involves addition of H to 
the end N of N 2 O which has a 10.3 kcal/mol barrier and a subsequent 1,3-hydogen 
shift to yield products via the NNOH saddle point which has a 16.3 kcal/mol barrier. 
Thus, the indirect pathway has a lower overall barrier and should be the more 
favorable process at lower temperatures. 
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The reaction of the 1 A excited state of NH with NO, reaction (5), has also been 
studied. The reaction of NH 1 A with H 2 has been studied by Fueno et al. [23]. These 
authors observed an initial approach geometry in which the N atom approaches 
perpendicular to the midpoint of the H 2 bond and the NH bond is perpendicular to 
the plane formed by the N and H 2 molecule. Fig. 2 shows the electronic structure 
for the lowest six surfaces of NH + NO for this geometric arrangement. Two of 
these surfaces arise from the 3 S“ state of NH and the remaining four surfaces arise 
from the 1 A state of NH. Of these surfaces the one which allows insertion is given by 
Eqn. (9). A limited study of this insertion process was made using a state averaged 
CASSCF calculation (averaged over the six lowest roots) varying only the N to NO 
bond midpoint distance ( R ) as shown in Fig. 3. From Fig. 3 curve crossings 
are evident between the diabatic states arising from NH 3 E _ state (repulsive) and 
the diabatic state arising from Eqn (9). If the H atom is allowed to dissociate 
from the initially formed HN-NO adduct, there axe four electrons of it' and four 
electrons of tv" symmetry with respect to the NNO plane. Thus, this structure does 
correlate with the ground state of N 2 0 and it is plausable that production of N 2 0 
via addition of NH 1 A to NO occurs via the third excited surface. 

V. Conclusions. 

The potential energy surface for NH + NO has been characterized using CASSCF 
gradient calculations to determine the stationary point geometries and frequencies 
followed by CASSCF/ internally contracted Cl calculations to determine the ener- 
getics. The CASSCF method is found to give reasonable geometries and frequencies 
at minima on the PES. At saddle points the geometries and frequencies are some- 
what less accurate as revealed by comparison of CASSCF and CCI results for the 
HN 2 and H0 2 PES’s. Studies of the HN 2 PES indicate that the CASSCF surface 
is distorted from the more accurate CCI surface mainly in the direction of the re- 
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action coordinate, while directions orthogonal to the reaction coordinate are better 
represented. 

In general the present results are in qualitative accord with the work of Melius 
and Binkley using the BAC-MP4 method. However, the present calculations have 
no empirical corrections while the BAC-MP4 method involves large corrections 
(~ 10 kcal/mol for single bonds and « 20 kcal/mol for multiple bonds). There 
are also some large discrepancies (as large as 8 kcal/mol) in detailed energetics 
between the present work and the BAC-MP4 results. 

Addition of NH to NO on a 2 A" surface involves barriers of 3.2 kcal/mol for trans 
orientation and 6.3 kcal/mol for cis orientation. The presence of these barriers is 
predicted from simple qualitative considerations and is supported experimentally 
by threshhold measurements for production of NH + NO in the reaction of hot H 
atoms formed by photolysis of an HI-NO complex in the experiments of Bohmer et 
al. 

Production of N 2 + OH products is predicted to occur via a 1,3-hydrogen shift 
from an initially formed cis HNNO adduct of planar 2 A' symmetry. This conclu- 
sion is supported by experiments of Petel-Misra and Dagdigian in which a marked 
preference for H(A') A doublet levels is observed in the OH product. 

Addition of NH to NO is shown to occur on an excited state surface and to 
involve a geometry in which the NH approaches the NO with the N perpendicular 
to the NO bond midpoint and with the NH bond perpendicular to the NO-N plane. 
The electronic structure of this complex correlates with H plus the ground state of 
N 2 O. This suggests that N 2 O observed as the primary product in the reaction of 
NH 1 A with NO by Yamasaki et al. occurs on this excited state surface. 
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Table I. Computed geometries and frequencies for fragments' 1 . 




calc. 

exp. 

NH 

?NH 

1.062 

1.038 


UJ 

3131 

3282.09 

NO 

*NO 

1.174 

1.1508 


(jO 

1889 

1904.3 

N 2 0 

r nn 

1.149 

1.1282 


TjVO 

1.204 

1.1842 



2267 

2223.7 


U?2 

1297 

1276.5 



568(2) 

589.2 

OH 

r oh 

0.988 

0.9706 


UJ 

3624 

3735.21 

n 2 

?NN 

1.123 

1.0977 


LO 

2315 

2358.027 


bond lengths in A and frequencies in cm 


-i 
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Table Ila. Comparison of CASSCF and CCI for HO^ 



H-0 2 


H0 2 



CAS/grad 

CCI 

CAS/grad 

CCI 

r oh 

3.55 

4.14 

1.86 

1.84 

r 00 

2.33 

2.29 

2.58 

2.52 

L H00 

116 

116 

102 

104 

oq 

596 i 

412 i 

3617 

3531 

(x >2 

1503 

1515 

1438 

1417 


501 

354 

1058 

1220 

E 

1.9 

0.4 

-39.9 

-51.2 


a Bond lengths in a 0 , angles in degrees, frequencies in cm and E in kcal/mol. 
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Table lib. Comparison of CASSCF and CCI for HNf 



h-n 2 


hn 2 



CAS/grad 

CCI 

CAS/grad 

CCI 

r NH 

2.49 

2.64 

2.05 

2.02 

*NN 

2.19 

2.17 

2.27 

2.25 

L HNN 

116 

123 

114 

115 


1788 i 

1387 i 

2577 

2577 

U> 2 

2004 

2027 

1735 

1931 

U?3 

892 

668 

1116 

1106 

E 

27.4 

15.8 

23.6 

5.6 


a Bond lengths in ao, angles in degrees, frequencies in cm and E in kcal/mol. 



Table Ilia. Computed frequencies and geometries ( 2 A' surface) 
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Table Illb. Computed frequencies and geometries ( 2 A" surface) 
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Table IV. Computed energies and zero-point corrections. 



Energy 

zero-point energy 

A E 

NH+NO 

-184. 81050(-. 85245) 

0.01143 

0.0 

H+NNO 

-184. 85597(-. 90221) 

0.01071 

-31.7 

cis HN-NO 

-184.78351 (-.84556) 

0.01457 

6.3 

trans HN-NO 

-184. 78893(-. 85093) 

0.01507 

3.2 

cis HNNO 

-184. 87106(-. 93448) 

0.01979 

-46.2 

trans HNNO 

-184. 87834(-. 94352) 

0.02036 

-51.5 

NNO-H 

-184. 80993(-. 87431) 

0.01142 

-13.7 

H-NNO 

-184. 82843(-. 88698) 

0.01189 

-21.4 

NNOH 

-184.81960(-. 87912) 

0.01350 

-15.4 

cis HNNO 2 A" 

-184. 83465(-. 89268) 

0.01746 

-21.5 

trans HNNO 2 A" 

-184. 83635(-. 89190) 

0.01732 

-21.1 
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Table V. Comparison of BAC-MP4 and CASSCF/CCI for NH + NO; Geometries and Frequencies 
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Table VI. Comparison of BAC-MP4 and CASSCF/CCI for NH + NO; Energetics 


Stationary Point 

BAC-MP4 a,i> 

CASSCF/CCI' 

H + N 2 0 

0.0 

0.0 

NNO-H 

19.8 

18.0 

H-NNO 

3.6 

10.3 

[NNOH] 

16.7 

16.3 

cis HNNO 2 A! 

-14.5 

-14.5 

trans HNNO 2 A' 

(-11.6) 

-19.8 


° relative energy in kcal/mol. 
b 6-31G** basis set 
c [4s3p2dlf/3s2pld] basis set 
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Figure Captions. 


Fig. 1. Schematic of the potential energy surface for NH + NO. The stationary point 
geomtries and vibrational frequencies are obtained at the CASSCF level, while the 
computed energetics are from Cl calculations ( See the text.). A region of crossing 
between 2 A' and 2 A" surfaces in the NH + NO entrance channel region is indicated. 

Fig. 2. The asymptotic electronic structure of the six potential energy surfaces 
arising from NO 2 H plus NH 3 S~ and NH 2 A. 

Fig. 3. Potential energy curves for the lowest six potential energy surfaces for NH 
+ NO as a function of R with the other geometrical parameters fixed (See the figure 
inset.). 
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